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RECONSTRUCTION  OF  PERIODIC  BOUNDARY  BETWEEN 

DIELECTRIC  MEDIA 


Jean  Chandezon1,  Anatoliy  Ye.  Poyedinchuk2,  and  Nataliya  P.  Yashina1,2 

1 Universite  B.  Pascal,  Clermont-Ferrand,  France 
2 IRE  NASU,  Kharkov  61085,  Ukraine.  Email:  nataliya@lin.com.ua 

Boundary  shape  analysis  and  reconstruction  is  important  area  of  research  [1  -4],  however  there  is 
a shortage  of  rigorous  approaches.  We  propose  a robust  and  clear  modification  of  the  known  C 
method  [5,6]  for  solving  the  wave  scattering  by  arbitrarily  shaped  surface.  This  approach  makes 
a reliable  basis  for  solution  of  recognition  problem:  the  reconstruction  of  surface  profile  and 
material  parameters  of  media  from  known  data  on  the  scattered  field. 

For  the  direct  problem  solution,  the  C-method  in  combination  with  a regularisation  has  been 
chosen  [5,6].  This  enables  us  to  reduce  the  original  2D  problem  of  linearly  polarized  plane  wave 
diffraction  by  an  arbitrary  boundary  of  dielectric  media  to  operator  equation 

aX  + AX  = C , (1) 


where  a is  regularizing  parameter  , A is  a self-adjoint  positively  defined  compact  operator  in 
corresponding  Hilbert  space  H = l2  x/2 . The  entries  of  unknown  vector  X are  coefficients  of 
scattered  field  that  is  expanded  over  eigenfunctions  of  C-method  [5,6],  Vector  C is  defined  by 
the  excitation  field.  For  boundary  shape  between  media  and  given  dielectric  parameters  the 
equation  (1)  can  be  solved  efficiently  by  means  of  truncation  method.  The  choice  of  regularising 
parameter  a may  be  done  on  the  base  of  generalised  residual  principle  [7,8],  The  input  data  for 

correspondent  inverse  problem  are  complex  amplitudes  R = {R„{^)X=-n  of  reflected 
propagating  wavesX  is  a wavelength.  We  suppose  that  this  data  are  known  in  certain 
wavelength  range  [A.;  ] . Besides  the  period  of  boundary  shape  and  dielectric  parameters  of 

media  are  also  known.  It  is  necessary  do  find  out  by  these  input  data  the  function,  describing  the 
boundary  of  two  media.  Let  a = ( an  be  Fourier  coefficients  of  this  boundary  function.  The 

solution  of  operator  equation  (1)  gives  the  mapping  that  associates  set  of  a-  with  set 

of  complex  amplitudes  R = (Rn(x))^  N ■ Thus,  on  certain  set  of  vectors  D,,  al2  the  non  linear 
operator 

F(<7,A.)=  7?(X.),  ^2) 

is  defined.  Having  found  out  from  (2)  the  Fourier  coefficients  a = (an  )“=_M , we  can,  by 


summation  of  Fourier  series  with  approximate  in  l2  -space  metric  coefficients  [7],  derive  the 
profile  function.  Formally,  the  scheme  can  be  outlined  as  follows.  Let  Y( l),  that  is  the  solution 
of  (2),  is  the  set  of  operator  F values.  Introduce  on  T(X)  the  norm  according  following 


formula  =^|/?(A,)|  cos®"  . Here  following  notations  are  used:  0„  are  angles  of 

diffracted  field,  0 is  angle  of  incident  field.  Consider  the  functional  that  is  given  in  definitional 
domain  D,,  of  operator  F 
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K-^Sj 


(3) 


=1  n=-P(k„) 


n=-0 


where  y > 0 , is  the  parameter  of  regularization,  R > 0 , p„  = ^K2m£,i±  -n2  , Km  = d/Xm  , d 

o , , 

is  a period  of  boundary,  ^e[^A2],  fl(v)  = • Vector  = \ayn )^=  , which 


n=-0 


provides  the  functional  (3)  with  minimum  we  consider  to  be  a solution  to  (2). 


Numerical  experiments.  On  the  basis  of  the  approach  developed  the  numerical  algorithms  for 
the  solving  (1)  and  (2)  have  been  implemented.  The  search  of  vector  a yp  is  organized  by  means 

of  regularized  quasi  Newton’s  method  with  step  adjustment,  using  only  fist  derivatives.  The 
minimum  residual  method  is  applied  for  the  choice  of  regularizing  parameter  y . Based  on  (1), 

we  simulated  input  data  Re (km ) = (fl* (Xm )fn=_N  ,m  = l,2...P  for  several  boundaries  between 
media  characterized  by  profiles: 

a,{y)  = h 
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periodically  continued  from  interval  [0,d]  onto  interval  (-co,+oo).  Parameters  d and  h 
satisfy  the  restriction  2nh/d  < 1 . The  wavelength  of  incident  plane  E polarized  wave  was 
varying  within  the  range0.5  < d/X  < 3.5.  Permittivity  of  the  first  medium  has  been  chosen  as 
s,  =1  and  of  the  second  one  as  s2  =2.25.  Permeability  of  both  media  is  p.  =1 . Functions 
a\  (v) , ai  (y)  are  chosen  for  they  belong  to  two  essentially  different  classes.  Namely,  function 
a2(y)  has  a finite  series  of  its  Fourier  coefficients.  In  the  contrary,  function  a^{y)  and  a3(y ) 

have  an  infinite  Fourier  series,  which  Fourier  coefficients  have  algebraic  type  of  decaying  only. 
Results  of  numerical  tests  are  presented  in  Fig.  1.  Solid  lines  correspond  to  the  exact  functions 

ai{y)d  -1>2,3  Lines  depicted  as  crosses  are  the  graphs  of  functions  and  aR2\y)  that 

have  been  defined  via  input  data  Re (Xm ) = (i?‘’ (km  ))^=_w  according  to  above  described 
algorithm.  As  they  almost  coincide  with  graphic  accuracy,  the  deviations 


10h-‘\a(y)-aRi(y)\  i = 1,2,3  are  presented  in  the  same  figures  as  doted  lines.  It  worth  to 


emphasize  that  maximum  absolute  value  of  deviation  essentially  decreases  with  value  of  points 
P increases  (we  remind  that  P is  a number  of  values  of  incident  plane  wave  wavelengths,  for 
which  the  input  data,  m=  1,2,.  ..  P,  have  been  calculated)  that  is  in  compliance  with  given  level 

of  noise  in  input  data  R(Xm)  = {R„{X„,))^=_n  . Rather  good  approximation  used  for  starting 

values  of  Fourier  coefficients  guarantees  the  shape  reconstruction  with  accuracy  10  2 -7(TJ . 
These  algorithms  are  efficient  tools  for  the  study  of  influence  of  input  data  errors  on  the 
accuracy  of  boundary  shape  reconstruction. 
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Figure  1.  Reconstruction  of  profiles  given  by  functions  a,  (y)  (a ),a2{y)  (b),  and  a3(y) 
(c)  for  the  values  P = 6,  e2  =2.25  ,0.5  < d/X  < 3.5,  2tc h/d  = 0.4. 


It  is  well  known  that  one  of  the  most  complicated  problems  in  solving  unstable  problems  of 
reconstruction  is  the  problem  of  matching  regularizing  parameter  a with  given  level  of  errors  in 
input  data  Remn  . We  have  demonstrated  that  such  tradeoff  can  be  obtained  be  means  of  residual 
method  [7,  8].  We  shall  demonstrate  this  statement  for  the  test  problem  for  reconstruction  shape 
of  the  surface,  described  by  the  function  a2  (y) . Basing  on  the  solution  to  direct  problem  we 

have  calculated  the  input  data  RLmn  = (A,„  )(l  + Y -Rcind)  for  various  levels  of  relative  error 


y . The  error  Rand  has  been  simulated  be  the  generator  of  random  numbers  with  normal 
distribution.  One  of  the  numerical  examples  is  presented  in  Fig.  2.  In  fragment  a)  you  can  see  the 
characteristic  behavior  of  relative  error  of  profile  reconstruction  has  been  estimated  according 

, , ;zk-<i2(> 

formula  5(a)  = | — : — | , calculated  for  numbers  of 


11+ Inf  If" 


a = 10~” ,n  — 1,2 ...N;  N < 10 . Here  aen  are  exact  values  of  function  a2(y)  Fourier 
coefficients.  As  it  is  clearly  seen  the  function  5(a)  has  pronounce  minimum,  that  appears  for  all 
considered  levels  of  error  y . In  pictures  b)  and  c)  the  results  of  boundary  shape  reconstruction  re 
presented.  Relying  on  our  numerical  experiment  we  can  conclude  that  reconstruction  can  be 
performed  with  best  accuracy  for  that  value  of  a , that  provides  minimum  to  function  5 (a). 

The  residual  method  can  provide  reliable  determination  of  optimal  value  of  a according  to 


relation  aopl  = sup(a  : Aa(y)<y).  Here  \ = 


iVm/7  ^ mr, 


V" 


II*. 

n 


is  a relative  residual  of 


J a 


input  data  R,im  and  i?“n,  that  are  the  results  of  the  solution  of  direct  problem,  calculated  for 


function  a'1  (y)  found  out  from  minimization  of  (3)  for  given  y . From  the  results  of  numerical 

experiments  we  can  see,  that  <\(y)  depends  on  a monotonically,  and,  thus  a , is  unique 

for  each  level  of  input  data  errors.  The  suggested  algorithm,  which  performs  reconstruction  of 
shape  of  periodic  boundary  between  two  dielectric  media  relying  on  information  about 
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diffraction  harmonics  that  are  known  within  certain  interval  of  wavelength,  requires  certain 
starting  approximation  for  Fourier  coefficients  of  function  a{y)  describing  boundary.  Such  an 
approximation  can  be  constructed  by  generalizing  results  of  [10]. 


Fig. 2 


point  b 


c) 
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